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NEUMANN PROBLEM FOR NON-DIVERGENCE 
ELLIPTIC AND PARABOLIC EQUATIONS WITH BMO^ 
COEFFICIENTS IN WEIGHTED SOBOLEV SPACES 


HONGJIE DONG, DOYOON KIM, AND HONG ZHANG 

Abstract. We prove the unique solvability in weighted Sobolev spaces 
of non-divergence form elliptic and parabolic equations on a half space 
with the homogeneous Neumann boundary condition. All the leading 
coefficients are assumed to be only measurable in the time variable and 
have small mean oscillations in the spatial variables. Our results can 
be applied to Neumann boundary value problems for stochastic partial 
differential equations with BMOa, coefficients. 


1. Introduction 


In this paper, we study Lp estimates for elliptic and parabolic equations 
in non-divergence form: 

aijDijU + biDiU + cu — Xu = f in 

—ut + aijDijU + biDiU + cu — Xu = f in (— 00 , T) x 
with the homogeneous Neumann boundary condition, where A is a nonneg¬ 
ative constant and R!^ is a half space defined by 

R+ = {x = (xi, • • • , Xrf) = (xi, x') : xi > 0}. 


We consider the equations in weighted Sobolev spaces with measures 
^d{dx) = dx and fi{dx dt) = x^~‘^ dx dt 

in the elliptic and parabolic cases, respectively, for some 6 G {d — 1, d—l+p)- 
Krylov [^] first studied Laplace’s equation and the heat equation in 
weighted Sobolev spaces Hjg and see Section 2 for precise defini¬ 

tions. After 1 ]^, there has been quite a few work on the solvability theory 
for elliptic and parabolic equations in weighted Sobolev spaces, for instance, 
1^, |l^, 1^|. In particular, the authors of |]K, K] studied second- 


see 


order parabolic equations with the Dirichlet boundary condition in weighted 
Sobolev spaces with leading coefficients having small mean oscillations. The 
motivation of such theory came from stochastic partial differential equations 
(SPDEs) and is well explained in |16]. 
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Recently, Dong and Kim Q studied both divergence and non-divergence 
type elliptic and parabolic equations on a half space in weighted Sobolev 
spaces with the Dirichlet boundary condition. The coefficients in Q are 
contained in a larger class than those in |l^. Namely, the leading coef¬ 
ficients are assumed to be only measurable in t and xi except an, which is 
measurable in either t or xi, where xi is the normal direction. Kozlov and 
Nazarov considered an oblique derivative problem for non-divergence 
type parabolic equations on a half space with coefficients discontinuous in 
t (and continuous in x) in a weighted Sobolev space. Their proof is based 
on a careful investigation of Green’s functions. In this paper, we extend the 
result in \14\ to a more general setting. Namely, the coefficients considered 
in this paper are measurable in the time variable and have small mean oscil¬ 
lations with respect to a weighted measure in the spatial variables. We call 
this class of coefficients BMOa,. The weight, for instance, for the parabolic 
case is n{dxdt) = x\~'^ dxdt, where 6 G {d — l,d — 1 + p). The condition 
0 G {d — l,d — 1 + p) is sharp even for the heat equation; see [^]. We 
note that the coefficients aij in Q also have small mean oscillations with 
respect to a weighted measure as functions of x' G (whereas, in this 

paper as functions of x G M'^), but the size of the modulus of regularity of 
is proportional to the distance to the boundary. See Assumption O] and 


Remark 2.2 


Since the counterexamples of Ural’ceva |^] and Nadirashvili [^|, particu¬ 
lar types of discontinuous coefficients have been considered for the solvability 
of equations. One type of discontinuous coefficients, which has been widely 
considered, is the class of vanishing mean oscillation (VMO) coefficients. 
The study of equations with VMO coefficients was initiated by Chiarenza, 
Frasca, and Longo |]|, ^. In []^| Krylov gave a unified approach to inves¬ 
tigating parabolic and elliptic equations in unweighted Sobolev spaces with 
coefficients that are measurable in the time variable and have small mean os¬ 
cillations with respect to the usual Lebesgue measure in the spatial variables 
(BMOa; with respect to the Lebesgue measure); see also |^. In fact, the 
coefficients in [^] are called VMOa, coefficients, but their mean oscillations 
in X do not have to vanish as the radii of cylinders go to zero. For more 
related work about Lp theory with BMOa, or partially BMOj, coefficients for 
parabolic systems and higher-order parabolic systems, we refer the reader 
to Biii and the references therein. 

Our proof is in the spirit of the approach by Krylov |^. The key point of 
such approach is to establish mean oscillation type estimates for equations 
with simple coefficients, i.e., coefficients are only measurable as functions of 
t. Then we apply a perturbation argument, which is well suited to the mean 
oscillation estimates, to deal with BMO^, coefficients. Finally we obtain 
the desired Lp estimates by applying the Fefferman-Stein theorem on sharp 
functions and the Hardy-Littlewood maximal function theorem in weighted 
Lp spaces. 
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Here one of the main steps is to get the mean oscillation estimates of D^u. 
For a simple equation 

-ut + aij{t)DijU = f 

in M X with the Neumann boundary condition Diu = 0 on {xi = 0}, 
we treat DDiu and D‘^,u separately. We estimate DDiu as follows. Differ¬ 
entiating the equation above with respect to xi, it is easily seen that Diu 
satisfies the divergence type parabolic equation 

-{Diu)t + Di{aijDj{Diu)) = Dif 

in R X with Diu = 0 on {xi = 0}. Therefore, we can apply a result in 
0 to obtain the mean oscillation estimates of DDiu. On the other hand, 
the estimates of are much involved. We treat the mean oscillations 

of D'^/U in the xi variable and x' variables differently. By integrating by 
parts and the Poincare inequality in weighted spaces, we manage to bound 
the mean oscillations of in the xi variable by the maximal functions 
of DDiu. For the mean oscillations in x' variables, we write the equation in 
the following form 

d 

-Ut + ^ UijDijU = f - + aij)DijU - auDuu, 

i=2 

which can be regarded as a non-divergence type parabolic equation in R x 
applying an interior estimate result without weights for 
where u is, as a function of x' E a solution of a non-divergence type 

equation in the whole space R'^”^ (see, for instance, |]T^]), we bound the 
mean oscillations of D^iU in the x' variables by the maximal functions of /, 
DDiu and D‘^,u. 

As an application of our results, in a forthcoming paper we are going to 
study non-divergence form SPDEs in weighted or unweighted Sobolev spaces 
with the Neumann boundary condition. A particular case is the solvability 
of SPDEs in the form 


du = {aijDijU -|- biDiU + cu — \u + f) dt + dw^ in (—oo, T) x 




where are independent one-dimensional Wiener processes, atj, bi, and 
c satisfy the same conditions as in the current paper, and /, Dg^, gk E 
Ihp,6»(—oo, r); see the definition of the 'hpfi space at the beginning of Sec¬ 
tion 2. We note that SPDEs in weighted Sobolev spaces with the Dirichlet 
boundary condition have been studied extensively in the past fifteen years. 
We refer the reader to |^, and the references therein. 

This paper is organized as follows. In the next section, we introduce 
some notation and state our main results. In Section 3, we obtain the mean 
oscillation estimates for D^,u and DDiu separately for a parabolic equation 
with simple coefficients. In Section 4, we prove our main theorem (Theorem 
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2. PRELIMINARIES AND MAIN RESULTS 


Throughout the paper we use, for example, the following Einstein sum¬ 
mation convention: aijDijU = Yhi j We introduce some notation 

used in the paper. As hinted in the introduction, a point in is denoted 
by X = (xi,--- ,Xrf), and also by x = (xi,x'), where x' G Recall 

= {x : xi > 0}. In the parabolic case a point in = R x is 

denoted by X = (t,x). Set = {(t,x) : xi > 0}. For r > 0, let 
be the open ball in R'^“^ of radius r with center x'. Denote 


Br{x) = Br{xi,x') = (xi — r,xi + r) x R'(x'), 
Qrit, x) = [t — ,t) X Br{x), 

R+(x) = Br{x) nR^J_, Q^{t,x) = Qr{t,x) 

For a G R, we use Qr{o) to denote 

Qr{Q, a, 0) = (—r^, 0) x {a — r,a + r) x 5^(0), 


and Qr = Qr(0). Similarly, we define Q^{a) and Q+. By a'^ we mean 
max{a, 0 }. 

Throughout the paper, we assume that the leading coefficients aij are 
bounded, measurable, and satisfy the ellipticity condition: 

aijiliflj > \aij\ < 1/(5 


for any rj G R'^, where <5 > 0 is a constant 

To introduce the function spaces used in this paper 
weighted Sobolev spaces 
integer 




introduced in jD]. 


we first recall the 
If 7 is a non-negative 


Hie = Hlei^i) = {u : x'^'d^u G Tp,0(R^) Va : 0 < |a| < 7}, 

where Lp^ 0 (R^J.)(= Lp^g) is a Lebesgue space with the measure Hd{dx) = 
x^~'^dx. For a general real number 7 , is defined as follows. Take and 
fix a nonnegative function C G (17^(0,00) such that 


00 

> 1 

n=—oo 


for all xi G R. For any 7 ,0 G R and p G (l,oo), let Hjg be the set of all 
functions u on R 1 such that 


|P = 


E 


e”'’||»(e”.)C(rci)||!;,p<oo, 


where || • is the norm in the Bessel potential space //^(R'^). For any 
a G R, let M“ be the operator of multiplication by (xi)" and M := M^. 
We write u G if M-^u G We set 

H;,(S,r) = Lp{{S,T),Hlg), Lp,0(S,r) = Lp{{S,T),Lp^g), 
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where —oo < S < T < oo. 

Our solution spaces are defined as follows. For the elliptic case, we set 

= {u:u, Du, D'^u G 

For the parabolic case, 


T) = {u. u, Du, D\ ut G Lp, 0 ( 5 , T)}, 
where —oo < S < T < oo. 

We also use the following Holder spaces. For a function f on D C 
define 


[f]a,b,V =■ sup 


f \f{t,x) - f{s,y)\ 
I |t — s|“ + |x — y\^ 


: {t,x), {s,y) G V, {t,x) 7 ^ {s,y 


where a,b G (0,1]. For a G (0,1], we set 


|,a,D = ll/l|oo,D + [f]^,a.V- 


The space corresponding to || • \\a/2,a x> is denoted by 
Throughout the paper, we use the weighted measures: 

lXd{dx) := {xi)^~'^dx, ix{dxdt) := {xi)^~'^ dx dt. 


where 9 G {d — l,d — 1 + p) and p G (1, 00 ). 

Now we state our regularity assumption on the leading coefficients. For 
a function g on denote 

b(T-)]B+(a;) = / 9{t,z) pd{dz)\pd{dy). 

^ ’ J B+{x) J B+{x) 

Then we define the mean oscillation of g in Q^(s,y) with respect to x as 

osCa;ig,Qtis,y)) = i b(T •)]!?+ 

J s—r"^ 

and denote 

g*= sup sup oscxig,Qt{s,y)). ( 2 . 1 ) 

Using the notation above with aij in place of g, we state the following 
regularity assumption on atj with a sufficiently small parameter p > 0 to be 
specified later. 


Assumption 2.1 (p). There exists a positive constant Rq such that 


A# := sup < p 


# 




Note that under this assumption, the coefficients a*-^ may not have any 
regularity with respect to t. 
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Remark 2.2. While we have a fixed size of the modulus of regularity Rq 
above, in Q the size of the modulus is proportional to the distance from the 
boundary to the location where the mean oscillations of are measured. To 
express this, one can replace R in (2.1) by yiR- This means, in particular, 
that the coefficients aij in |Q] are allowed to be much rougher near the 
boundary than those in this paper. 


For lower-order terms, we assume that the coefficients bi and c are only 
measurable (without any regularity assumptions) and bounded so that 

\bi\, |c| < K 

for some constant K > 0. 

The following theorems are our main results, the first of which is the 
unique solvability of parabolic equations. 

Theorem 2.3. Let T G (—oo, oo], 1 < p < oo, and 6 G {d — l,d — 1 + p) 
be constants. Then there exist constants p = p{d,6,p,9) > 0 and Aq = 
Xo{d,6,p,9, K, Rq) > 0 such that under Assumption (p) the following 
assertions hold. 

(i) Suppose that u G (—oo,T) satisfies 

— ut + OijDijU + biDiU + cu — Xu = f (2-2) 

in (—oo,r) X with the Neumann boundary condition Diu = 0 on {xi = 
0}, where f G hp^ 0 {—oo,T). Then we have 

||^t||p,6» +11-^ ^Diu\\pfi + \\D‘^u\\p^g + y/X\\Du\\p ^0 +X\\u\\p ^0 < CoII/Up,6* (2.3) 

provided that X > Aq, where || • = || • ||Lpe(-oo,T) Cq = Co{d,6,p,9). 

In particular, when bi = c = 0 and aij = aij(t), we can take Aq = 0. 

(a) For any A > Aq and f G Lp^p(—oo,T), there is a unique solution 
u G VFpg(—oo,r) to the equation (p^). 

We now present our results for elliptic equations, where the coefficients are 
independent of t. Since Assumption |2.1| (p) does not concern the regularity 
of coefficients in t, we still require the coefficients to satisfy Assumption |2.1| 
(p)- 

By adapting, for example, the proof of |l^. Theorem 2.6] to the results 
above for parabolic equations, i.e., by regarding elliptic equations as steady 
state parabolic equations, we obtain the following theorem for elliptic equa¬ 
tions. 


Theorem 2.4. Let 1 < p < oo and 9 G {d—l,d—l+p) be constants. Then 
there exist constants p = p{d,6,p,9) > 0 and Aq = Xo{d,6,p,9, K, Rq) > 0 
such that under Assumption (p) the following assertions hold. 

(i) Suppose that u G W‘^q{W\_) satisfies 

OijDijU + biDiU + cu — Xu = f 


(2.4) 
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in with the Neumann boundary condition Diu = 0 on {xi = 0}, where 
f £ Then we have 


||M ^Diu\\pfi + ||D^u||p^6i + \/A||-Dm||p^6i + A||tt||p^e < CoH/Hp^e 


provided that A > Aq, where || • ||p^6» = || • llLpg(K^) C'o = Co{d,5,p,6). In 
particular, when 6^ = c = 0 and aij are constant, we can take Aq = 0. 

(a) For any A > Aq and f G Lp^ 0 {M.'^), there is a unique solution u G 
to the equation (|2.4|) . 

3. EQUATIONS WITH COEFFICIENTS INDEPENDENT OF X 
In this section, we deal with equations in the form 



ut - aij{t)DijU = f in 
Diu = 0 on {xi = 0}. 


(3.1) 

(3.2) 


Note that now the coefficients aij depend only on t. Let us state several 
technical lemmas. The first one is Hardy’s inequality, which can be found 



Lemma 3.1. Let p G (1, oo), 9 G {d—l,d—l +p), and v G such 

that V = 0 on {xi = 0}. Then 



for any r G (0,oo], where C = C{d,p,9). 



ding theorem. 

Lemma 3.2. Let p G (l,oo) and a G (0,1) be constants. Assume that 

u G satisfies (|3.lD in with f = 0, and u = 0 on {xi = 0}. 

Then we have 

\\Du\\ca/2,a(^Q+'^ < C\W\\ Lj,{Q+)' 

where C = C{d, 6 ,p,a). If Qi and Q 2 replace Qf and Q\, respectively, then 
\\Du\\cc/2,^(q^^ < C||Dn||i^(Q2)- 

For a domain T) C we denote (n)© to be the average of n in P with 

respect to the measure ^{dxdt) = Xi~‘^ dxdt. Precisely, 



where 


Using Lemma |3.2|, we obtain the following mean oscillation type estimate. 







H. DONG, D. KIM, AND H. ZHANG 


Corollary 3.3. Letp E (1, oo), a E (0,1), yi E [0,1], and 6 E {d—1, d—l+p) 
be constants. Assume that u E satisfies 

ut - aij{t)DijU = 0 in 


with the Dirichlet boundary condition u = 0 on {x\ = 0}. Then there exists 
a constant C = C{d, 6,p, 9, a) such that for any r < 1, 


{\Du - < Cr^P [\Diu\P pidxdt). 

Proof. Prom Lemma |3.2| , we obtain 

Since 9 — d — p < 0 and by Lemma |3.1| , we have 

^P < (7 / dx dt < (7 f \Diu\P p{dx dt). 

jQi JQ.t 


\u\ 


LpiQt) 


Combining the two inequalities above, we prove the corollary. 


□ 


Before we state the next theorem, we introduce a function space. For 
A > 0, we denote u E oo,T) if 

'/Xu, M~^u, L>u E Lp^ 5 i(—oo, T), 

and Ut E oo, T) + ^/X]Lp^ 0 {—oo,T). Let us write 'Hpg{—oo,T) if 

A = 0. By Remark 5.3], one can find f,g = {gi,... ,gd) E Lp^ 0 (—oo,T) 
such that 

Ut = Digt + y/Xf 

in (— oo,r) X We set 

II^II^iA(_oo,t) = {^lkllp,0 + \W~^u\\p^e + ||L»u||p,e 

+ lbllp ,0 + ll/llp, 6 » : Ut = Digt + \/A/|. 

Now we state a special case (a*-^ = aX^{t)) of [§, Theorem 3.9], where a*-^ 
are allowed to be merely measurable in (t,xi) except that or 

= a^^(xi). Note that in the theorem below there is no specihcation of 
the boundary condition, but functions in the solution space 'Hp’g(—oo,r) 
necessarily satisfies u = 0 on the boundary. Hence the theorem is about the 
Dirichlet boundary value problem for divergence type equations in weighted 
Sobolev spaces. 

Theorem 3.4. Let T E (— oo, oo], A > 0, p E (1, oo), 0 E (d — 1, d — 1 +p), 
and u E oo,T) satisfy 

Ut - DfioijDju) + Xu = Digt + / 


(3.3) 
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in (-oo,r) X where g = (91,92, ■ ■ ■,gd), g,f & hp^g{-oo,T), and / = 0 
if X = 0. Then we have 

uWpfi + \\M ^u\\p^0 + ||-Dn||p^ 6 ) < Cll^llp^e + -^WfWp^g, 

where C = C{d,6,p,6). Moreover, for any g,fG hp^g(—oo,T) such that 
f = 0 if X = 0, there exists a unique solution u G 'Hpg{—oo,T) to the 
equation (|3.3D . 

Next we consider (H-Q with non-vanishing right-hand side. 



Lemma 3.5. Let p G (1, 00 ), 0 G (d—l,d—l+p), a G (0,1), r > 0, k > 32, 
and yi > 0. Assume that f G lLp^g(Qfj.(yi)) and u G is a solution 

of 


Uf a-ij D{jU — f 

in Q'f^{yi) with the Neumann boundary condition Diu = 0 on {xi = 0}. 
Then we have 


(|CD,u-(DAa)Q+(„.)l%+(,.) 

where C = C(d,5,p,6,a). 

Proof. Denote v = Diu. Then v satisfies 

vt - Di(aijDjv) = Dif 

in with the Dirichlet boundary condition n = 0 on {xi = 0}. By a 

scaling argument, it is sufficient to set nr = 8 . We consider two cases. 

Case 1: yi G [0,1]. Since r = 8 /k < 1/4, we have 

Qt(yi) c Q+ c g+ c Q+(yi). 

Let 7 ? be a smooth function with support in (— x B,^r(yi,^) 
and y = 1 in ^ 4 . By Theorem |3.4| , there exists a unique solution w G 
oo, 0 ) = Ti.^ q(—oo,C) to the equation 

wt - Di(aijDjw) = Di(fT]) 

in (—00,0) X satisfying 

\\Dw\\p^e < C\\fr]\\p^e- 

Due to the definition of 9 , this implies 

\\Dw\\%<Cp{Qtr(yi)){\f\%i^^y,y (3.4) 

By a standard mollification argument (see, for instance, [^, Theorem 4.7]), 
we may assume that w is smooth. Let u = v — w, which is also smooth, and 
satisfies ii = 0 on {xi = 0 } and 

Ui Di(aijDj'uj — 0 
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in Q 4 . By Corollary |^, we have 

(l^«- ^ [ ^\D,u\Pnidxdt). (3.5) 

J Q 4 

Combining (|3.4|), (^), and the triangle inequality, we reach 
{\DDlU - (^^l«)Q+(j/i)r)Q+(yi) 

< C{\Du-{Du)Qtiy,)\nQny.) + ^^\^^\^hn^ 

< Cr“?’(|Z)iir)g+ +C/r(Q+(yi))-V(Q;i'.(yi))(l/r)Q+(,o 

< Cr“?’(|DDinr)Q+ +C^(Q+(yi))-V(Q+(yi))(l/r)Q+(,,)• 

Bearing in mind that /r((5^(yi)) > where a = max{(i, 0} and r = 

8 /k, we prove the lemma for Case 1. 

Case 2: yi > 1. This is essentially an interior case. Since r = 8 /k < 1/4, 
we have 


Qtiyi) = Qrivi) c Qi/4(yi) c Qi/ 2 (yi) c Qtriyi)- 

As in Case 1, we take a smooth function r] with support in (—x 
BKr{yi,0) and 77 = 1 on Qi/ 2 {yi)- Then by Theorem |3.4 there exists a 
unique solution w G 'Hp 0 (—oo,0) of the equation 

wt - Di{aijDjw) = Di{fT]) 
in (—00,0) X M^, satisfying 

\\Dw\\p^e < C\\fr]\\p^g. 


Then we get 

\\DM\li < Cf<(<3t(!/l)) (l/l’’W,,„) ■ (3.6) 

For the same reason as before, we may assume that w is smooth and so 
\s u = Diu — w. It is easily seen that 


ut-aijDijU = 0 in Qi/ 2 (?/i)- 
By Lemma p.2| , we have 

\\Du\\c^/^ 

From this and ( |3.6|) , 

{\DDlU - {DDlU)Q^(^yj^)\P)Q^(y^) 

< Ci\DU - + CilDwinQ^t^y,) 

< + Cy{Qriyi))-^y{QUyimf\%Uy.) 

< + c^(Q.(2/i))-V(Q+(yi))(l/r)Qt(.0- (3.7) 

Since for any xi G ( 7/1 — 1/2, yi + 1/2), 
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where Ci ^2 
that 


Ci,2{d,9), by (^), (^), and the triangle inequality, it holds 


{\DDiu - 

< Cr^P{\DDMnQ^,,iy,)+C^i{Qr{yl))-^^i{Qtr{ym^ 

Taking into account that ^i.{Qr{yi)) > N^{Q'^{vi))^ we prove the lemma 
for the second case. □ 


It remains to estimate DijU with i,j > 1. Let us first state a Poincare 
inequality in weighted Lp spaces. 

Lemma 3.6. Let a E M+, a E (—l,oo), r E (0,oo), p E [l,oo), ^{dx) = 
{xi)°‘dx, X E and u E Then 


f f 

JB+(a) Je 


Br (a) J Br (a) 


\u{x) — u{y)\^ v{dx) v{dy) 


< 


Cr^i'{B^{a)) [ \Du\^ v{dx) 

JB+ia) 


(3.8) 


where C = C{d,p,a). 


Proof. When a. E [0, oo), the inequality is proved in |13, Lemma 4.1] with 
a missing constant depending on d. For the sake of completeness, we here 
present a proof when a E (—l,oo). Since the weight is with respect to xi, 
we only prove (^) for d = 1. In fact, to prove (|3.8D for d > 1 we just need 
to combine the case when d = 1 with the unweighted Poincare inequality. 

Due to scaling, it suffices to prove ( ^.^ ) with r = 1. We further assume 
a E (0,2). Indeed, if a > 2, the inequality (|3.8D is equivalent to the usual 
Poincare inequality without weights. For each x,y E ((a — l)'’“,a + l), by 
Holder’s inequality. 


rmax{x,yj p 

u'{z) dz 

' mm{x,y} 

pjq ^max{x,y} 

< c 


r 

\u{x)-u{y)\P= / 

J mi 


min{iE,y} 
CL -\-1 


u'{z)\Pz^ dz 


/ \ p/q J- 

< C (y^-^^/P + x^-^^/P] / \u'{z)\Pz^ dz, 

^ ’ d(a-l) + 

where 1/p + l/g = 1 and C = C{p, a). Then, to conclude (|3.8|) , we integrate 
the above inequalities with respect to u(dx) and ^{dy), and use the fact that 
a E (0,2). □ 


In the sequel, we denote the standard parabolic cylinder in as 
Qi+\X) = {t-r\t)xBr{x), 
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where Br{x) is the Euclidean ball in with radius r and center x. For a 
function / on we dehne the average of / in without weight 

as 

{f)n‘i+ltv\= -'^11 - / f(t,x)dxdt. 

Theorem 3.7. Let p G (l,oo) and u G satisfy 

—ut + aij{t)DijU = f in 

Then there exists a eonstant C = C{d,5,p) such that for any k > A, r > 0, 
we have 

Proof. See |^, Theorem 5.1]. □ 

To estimate D^/U, we introduce the following notation. For X)i > 0, set 
= ((t)i - r)+, t)i + r) C M, pi{dxi) = x\~'^ dxi, and 

r'c)i+r 


Lemma 3.8. Letp G (1, oo), 9 G (d— 1, d— l+p),r > 0 ,k> 32, and t)i > 0. 
Assume that f G 'hpfi{Qfj.{X)i)) and u G is a solution of 


Ut CLijDijU — f (3.9) 

in Qtri^i) dliu = 0 on {xi = 0}. Then there exist constants C = 
C{d,6,p,6) and a = a{d,p,6) > 0 such that 


i.\rl.u (r'^'*i)Q+(,,)r)Q+(,|,) 

< Ck-u\dU%u,.) + + (I"Oi“I")qJ.(,.))' 


Proof. By scaling, we may assume that nr = 8. In this case r < 1/4 because 
K > 32. Let r] = ri{t,x') G (^“(M x with a unit integral such that 

supp(r/) C Qf and 

\Dx'r]\ <C{d)r~'^~‘^ on Qf, (3.10) 

where 

Qf := (—r^,0) X B^~^ := (—r^,0) x {|x'| < r}. 

We consider two cases. 

Case (i): t)i G [0,1]. By Holder’s inequality and the triangle inequality, 
it is easily seen that 


f \Dl,u-{Dl,u)g+,.\P p{dxdt) 

J Qqn (t)l ) 


< c 


+ f \D‘^,u{t,x) — D‘^,u{s,y)\P p{dx dt) p{dy ds). (3.11) 

4Q+(oi) 4Q+(t,i) 
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Since, by the triangle inequality, 
\Dl,u{t,x) - Dl,u{s,y)\P 


< C 


+ C 


+ C 


{t,x) — / ri{a, z')D‘^,u{(j,xi, z') dz'da 

JQi 

is,y)- / r]{a,z')Dl,u{a,yi,z)dz' da 

jQi 

/ r]{a, z'){Dl,u{a, xi,z') — D‘^fu{a, yi,z')) dz da 

jQi 


DIiU 


the right-hand side of ( |3.11D is bounded by C{I + 11), where 
I : = 


II : = 


^ D'^,u{t,xi,x') — / r]{s,y')D‘^,u{s,xi,y') dy' ds y{dxdt), 

JQt+i) JQi 

f f / {Dl,u{t,xi,z') — Dl,u{t,yi,z'))rj{t,z') dz' dt 

J Bl+{x)i) J Bl+(<oi) JQi 

■ y.i{dxi) midyi). 

Let us now estimate I and II separately and first consider II. By inte¬ 
grating by parts and Holder’s inequality, 


iQi 


xi, z') — D^iu{t, yi, z'))r]{t, z')dz' dt 


'Qi 


xi,z') - D^/u(t, yi,z'))D^:r]{t, z')dz' dt 


< / \D^iu{t,xi, z') — D,^iu{t,yi, z')\^dz' dt ■ { \D,^ir]\'^dz'dt)<i, 

JQi JQi 

where 1/p -|- 1/q = 1. Prom ( |3.10| ), we have 

\D^,r]\'idadz'y < Cr-^'^+^+P). 


Qi 

We plug the two inequalities above into II to achieve 

II < 

-/ / / \Dx'u{t,xi, z') — Dx>u{t,yi, z')\^dz' dt 

J Bl++i) J Bl+{x)Q JQi 

■ iJLi{dxi) yLi{dyi). 

Applying Lemma with d = 1, we get 

II <C f \DiD^,u\P y{dx dt) < CK'^^^+^{\DDiu\P)n+ ) 

JQt+i) 

because |Q/'(r)i)| > 

Next, we estimate I, which can be written as 


I = 


+ / ri{s,y'){D^,u{t,xi,x) — D^,u{s,xi,y')) dy' ds fi{dxdt). 

JQt+i) JQi 
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By Holder’s inequality, 
I < 


+ / \D‘^,u{t,xi,x') — D‘^iu{s,xi,y')\^ dy' ds iJ,{dx dt) 

jQti'oi) JQi 

P 

( [ \r]\'^dy'ds\ . 




Since 


we have 
I <C 


yj \y\^dy'dsj <Cr-^'^+^\ 

7 - + \Dl.,u{t,xi,x') — D'^,u{s,xi,y')\^ dy' ds yL{dx dt) 

7q+(15i) JQd 

C / I{xi) x\~'^ dxi^ 

J Bl+{r,A 


(3.12) 


where 


I{xi) = P P \D‘^iu{t,xi,x') — D‘^iu{s,xi,y')\^ dy' ds dx' dt. 

Jq$ Jq$ 


We now estimate I{xi) by writing the equation (^) as 


Ut + ^ ^ CLijDij 

*d>i 


u = — 


auDiiu - + aji)DijU + /. 

j=2 


Here, for each fixed xi E (0,oo), we regard u{t,xi,x') as a solution to the 
above equation defined in M x Then thanks to Theorem ^ with d in 

place of d + 1 and the triangle inequality, we get 


I{xi)<C + \D'^,u{t,xi,x') — {D‘^,u)Qd{xi)\^ dx' dt 

JQI: 


< CK-n\Dl,nr)g,Jx,) + CK"‘{|/|f)e.„(rCi) + CK.‘‘+H\DDM’’)QiMl), 

where, for fixed xi, (xi) is the unweighted average of g with respect 

to (t, x') in the d dimensional parabolic cylinder with radius nr. We plug 
the inequality above into ( 3.12D to obtain 


/ < 


J Br+im) \ 

+ K'^+^(|T>Zliu|^’)Qd^(xi)^x?“^dxi. 


(3.13) 
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Bearing in mind that where = maxjd, 0}, and 

r = 8 /K, we obtain from (|3.13D that 

By the definition ofi^,p + d — C — 1 >0. Combining the estimates of I and 
II, we complete the proof of Case (i) with a=p + d — C — 

Case (ii): t)i > 1. Set v = D^u, k = 1 ,... ,d, and note that v satisfies 
the divergence type equation 

vt - Di{aijDjv) = Dkf 

in i?Kr(t)i). Since this is an interior estimate, we do not care about the 
boundary value of v on {xi = 0}. Then we repeat the second part of the 
proof of Lemma with Dk in place of Di. The lemma is proved. □ 


4. 


PROOF OF THEOREM 


2.3 


In this section, we deal with operators with coefficients depending on 
both X and t. We denote L = aijDij and assume p > 1, A > 0, and 
0 G {d — l,d — 1 + p). 

First we need the following lemma. 


Lemma 4.1. Let 0 < r < R < oo and X = {t,x), Y = {s,y) G Then 

p{Qt{X)) <C{d, 9) p{Q+{Y)) (4.1) 

provided that Q^{X) n Q^{Y) / 0. 

Proof. By scaling, without loss of generality we may assume that R = 1. 
We then consider two cases depending on yi. 

Case 1: yi > 10. In this case, we have 

xi > 8, 1 < max{yi + 1, xi + r}/min{yi — 1, xi — r} < 11/7, 

which allows us to use the Lebesgue measure in comparing p {Qf{X)) and 
p (Q)''(T)). Therefore, (^3) clearly holds since r < 1. 

Case 2: yi G [0,10]. In this case, we have xi G [0,12). It is then easily 
seen that p{Qf{Y)) is bounded from below by a constant c{d,6) > 0 and 
p{Qf{X)) < p{Qf{X)) is bounded from above by a constant C{d,6) > 0. 
Thus, (|^) still holds. 

The lemma is proved. □ 


The following two lemmas are mean oscillation estimates for the operator 
aij{t,x)Dij. We prove them by using the mean oscillation estimates for 
aij{t)Dij proved in Section combined with a perturbation argument. 

Lemma 4.2. Let R > 0, k > 32, and /3, /3' G (1, oo) satisfying 1//3 + 1//3' = 
1. Suppose that u G is compactly supported on Q^{Xi), where 
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Xi = {i,x) G Moreover, Diu = 0 on {xi = 0} and f := —Ut + Lu. 

Then for any r > 0 and Y = (s, y) G we have 

(|DOi„- (DO,u)g+|j./)g+,y, < Co^-f/^dDDiuing+^ir) 

+ (4.2) 

where Cq = C'o(d, 5,p, 9) and Ci = Ci{d, 6,p, 9, /3). 


Proof. Throughout the proof, we assume Q^{Xi) r\Qf^{Y) 0. Otherwise, 
(4.2) holds trivially. Fix a 2 ; G and set 

LzU = aij{t, z)Diju{t, x). 


Then we have 


-Ut + LzU = f + {aij{t,z) - aij{t,x))Diju{t,x) := f. 


It follows from Lemma 3.5 with a = 1/2 and a translation of the coordinates 
that 


{\DDiu {DDiu)q+^y)\^)q+{y) 

< C„K-f''"(|DO,tt|l’)o+,(y, + C„K""+"(|/|'’)g+ (V), (4.3) 


where Cq = Co{d, 5,p, 9). By the definition of /, the triangle inequality, and 
the fact that u vanishes outside Q~^{Xi), we have 




(4.4) 


where 

h = — 0.ij{t,x))XQ+(^Xi)^L'^\^) Q+,.{Y)^ 

and Xg+(Xi) indicator function of Q^(Xi). Denote B~^ to be Bf^{y) 

if Kr < R, or to be i?/^(x) otherwise. Define in the same fashion. We 
note that 

p{Q+)<C{d,9)p{QUY)) (4.5) 


provided that Q~^{Xi) n Qfj.{Y) / 0. It is obvious if nr < R, i.e., Q~^ = 
Qf^{Y). If Kr > R, the inequality is proved in Lemma 4.1 . 

Combining (O) and ( [4.41) and taking the average of each term with re¬ 
spect to 2 ; in B'^, we reach 


{\DDtu-{DDtu)Q+^y/)Q+^y^ < CoK-P/\\DDtu\nQt^^y^ 

+ CoK‘^+<^+\mQt^^Y)+CoK^^'^^ S IzPdidz). (4.6) 

J 
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Since u vanishes outside by Holder’s inequality, we get 

/ IzHd{dz) 

J B+ 


triY)) / b + Iq 


KQ 

■ fi{dx dt) fid{dz) 


QMY)nQ+{Xi) 


\{aij{t,z) - aij{t,x))Diju\P 


< 




1 /* /* 

FTTHT / (/ \aij{t,z) - aij{t,x)fP ^i{dxdt)] Hd{dz) 

rr(b)) JB+ ^Jq+ ^ 


Md 


(4.7) 


/Qt(H)nQ+(Vi) ’ 

By the boundedness of aij, Holder’s inequality, and the definition of osca;, 
we have 


/ + (/+ - aij{t,x)f'P ^{dxdtfj Hd{dz) 

<Ci(^J J \aij{t,z) - aij{t,x)\id{dxdt) fid{dz)^ 

< Cl (/x((5+)osca;(aij,Q+))^^^ , 

where Ci = Ci{d,6,p, /3). From ( |47D , ( |4)^ , and (^), we obtain 


/ hpd{dz)<Ci{A*)YP'{\D‘^uf^ 
J B+ 


'p \^/0 

QtriYY 


Combining (|4.6|) and (4.9), we get (4.2). The lemma is proved. 


(4.8) 

(4.9) 
□ 


Following exactly the proof of Lemma ^ with Lemma |3.8| in place of 
Lemma we obtain the lemma below. 


Lemma 4.3. Let i? > 0, a > 32, and 13, j3' € (1, oo) satisfying If(3 + If 13' = 

1. Let u E C')^(M)^'^^) he compactly supported on Q^{Xi), where Xi E 
Moreover, Diu = 0 on {xi = 0} and f := —ut + Lu. Then for any r > 0 

and Y = (s,y) E we have 

{\D„.,u — {D„.,u)q+^y)Y’)q+{y) 

< Ct,^-“(\Dl,nV)Q*,^y, + C„K‘‘+'’+2(|/r + |ODia|'')gi„|l') 

where Cq = Co{d, 6,p, 9), Ci = Ci{d, 6,p, 9, /3), and a = a{d, 9,p) > 0. 


Next we recall the Hardy-Littlewood maximal function theorem and the 
Fefferman-Stein theorem on sharp functions. Let 

Q = {Q+{z) :z = (t,x) EM^,r >0}. 
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For a function g defined on the weighted (parabolic) maximal and 

sharp functions of g are given by 

Mg{t,x)= sup f \g{s,y)\g{dyds), 

QeQ,{t,x)eQ J Q 

g*(t,x)= sup f \g{s,y) - {g)Q\fj,{dyds). 

QeQ,(t,x)£Q J Q 

For any 6 > d — 1 and g G we have 

ll■^fl'llLp,s(R^+l) — ^ll5'llLp,e(R^(+^) ’ 

where p G (1, oo) and C = C{d,p, 9). The first inequality above is known as 
the Fefferman-Stein theorem on sharp functions and the second one is the 
Hardy-Littlewood maximal function theorem, for instance see ]l^, Chapter 

3] 

Now we are ready to prove our main theorem. 


Proof of Theorem By the method of continuity, it is enough to prove the 
a priori estimate (|2.3D . Moreover, since the set of functions in C°°{{—qo, T) x 
vanishing for large (t, x) is dense in (— oo, T), we only need to prove (|2.3|) 

for infinitely differentiable fnnctions with compact snpport. In this case, the 
proof of (^) can be divided into several steps. 

Step 1: We consider A = 0, = c = 0, T = oo, and u G van¬ 


ishing outside for some Xi G where Rq is from Assumption 

|2.1[ Let K > 32 be a constant to be determined later. We fix g G (l,p) and 
(3 G (1, oo), depending only on p and 6 , snch that fdq < p and 9 < d — 1 + q. 
Let 13' be such that 1/[3 + 1/[3' = 1. By applying Lemma 4.2 with q in place 
of p and using Assumption 1^, we obtain 


{DDiu)*{Y) + CoN('^+'^+2)/57W^/'?(|/|'')(y) 


for any Y G where Cq = CQ{d,5,p,9). This estimate, together with 

Fefferman-Stein theorem on sharp functions and the Hardy-Littlewood the¬ 
orem on maximal functions, gives 

\\DDiu\\p^0 < C'||(ZIiAirt)^||p^6) 

< CoK-^^^DDiu\\p,e + CoK^''+^+^^/^f\\p,e + p^/^^''i^\D\\\p,e. 

(4.10) 
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In the same way, we apply Lemma 4.3 instead of Lemma 4.2 to obtain the 
estimate of D'^,u: 

\\Dlu\\p,e < CoK-^/nDl,u\\p,e + + \\DDiu\\p^e) 

+ (4.11) 

By choosing n sufficiently large, from ( |4.10 ) and (4.11) we obtain 

\\DDiu\\p^e < Co||/||p,0 + Cop^/(^'^)||Z?2^(||p^,, (4.12) 

\\Dl,u\\p^e<Co{\\f\\pfi + \\DDiu\\p,e) + Cop^^^^'‘'H^^ (4.13) 

We combine ( |4.12D and ([4.13| ) together to get 

WDMpfi < + 

Taking p sufficiently small depending only on d, 6,p,9, we arrive at 

\\DMp,e < Co\\f\\p,e. (4.14) 

On the other hand, since 

Ut — Qij Dij U /, 

we have 

||'Wt||p,6» < C'oll/llp.e + CoWD'^uWpfi. (4-15) 

By Hardy’s inequality (Lemma P) with r = oo, 

\\M~^Diu\\pfi < CQ\\Diiu\\pfi. (4-16) 

Combining ( 4.14| ), (|4.15|) , and ([4.16| ), we obtain ( |^) for A = 0. 

Step 2: We remove the assumption that u is compactly supported in 
(5^^(Xi). By a standard partition of the unity argument with ( 4.14|) ~ ([4.16| ) 
(cf. Ijl^, Theorem 1.6.4]) we see that 

||^t||p,0 + 11-^ + ||iA^u||p^0 < Coll/llp^e + Ciljujlp^e, (4-17) 

where Cq = Co{d,5,p,6) and Ci = Ci{d,5,p,9, Rq). 

Step 3: We still assume that bi = c = 0 and T = oo, but A is not 
necessarily zero. In this case, we follow an idea of S. Agmon. Since 

-ut + aijDijU = f - Xu, 

by ( 4.17 ) we have 

+ 11-^ ^Diu\\p^g + ||iA^u||p^0 < Coll/ — Att||p^0 + Cillttllp^e 

<Co\\f\\p,9 + {CoX + Ci)\\u\\p,e. 

Hence it is sufficient to show that for large A 

A||u||p,0 < C'o||/||p,6). 

We pick a function Ciu) £ C'(/’(M),C ^ 0 and introduce the following 
notation 

z = {x,y), u{t,z) = u{t,x)({y)cos{VXy), Lu = L{t,x)u{t, z) + Uyy{t, z). 
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Finally, set = {{t, z) : xi > 0} and 

Mr{z) = (xi — r, Xi + r) X Br{x', y), Qr{t, z) = {t — r'^,t) ^ 

B+ (z) = Br( 2 ;) n {xi > 0}, Q+ (t, z) = {t — r'^ ,t) x i^)- 


For any r > 0, Z = (t, z) G and a{t), we have 

f \aij{t,x)-a{t)\fj,d{dx)dy <C f \aij{t,x) - d{t)\y,d{dx), 
Jm+{z) J B+{x) 

where C = C{d). In particular, by the definition of the osc^; and setting 

d{t) = f aij{t,x) yLd{dx), 

J B+{x) 

we have osCz{aij,Q^{t, z)) < CosCx{aij,Q^{t,x)). By a simple calculation, 
Lu = f cos{VXy)C{y) + u(" cos{'/Xy) — 2VXuC' sm{VXy) := /. 


By (4.1?) with u instead of u in dimension d + 2, 


\D 




(4.18) 


where 


Wpfi 


is the weighted Lp norm in with respect to y,{dxdt) dy. 


Note that for A > 1, 

Cl < / \C{y) cos{VXy)\P dy := A < C2, 


(4.19) 


where ci and C2 are positive constants independent of A. Moreover, A is 
bounded from above with ( replaced by any derivatives of C and cos(\/Ay) 
replaced by sin(\/Ay). Then from ( 4.18 ), we have 

lll-^ra'“lllp ,0 — lll-^^^lllp,e — C'o||/||p,6» + Ci{l + A/A)||u||p,e- (4.20) 

Since 

DyyU = u{t, x)C,"{y) cos{V^y) — 2'/Xu{t, x)C'{y) sin(\/Ay) 

- Xu{t, x)C{y) cos{VXy), (4.21) 

Combining (|4.19D , (|4.20|) , and ([4.2l|) , we have 

< C'oll/llp,^ + Ci{l + \/A)||m||p^6i. (4.22) 

After choosing A sufficiently large depending on d,6,p,9, Rq to absorb the 
term of u on the right-hand side of (4.22|) to the left-hand side, we have 


All-Up,e < Co\\f\\p,e, 

where Cq = Co{d,S,p,9). By the interpolation inequality 

^/A||-^“IIp,6» < Co\\D'^u\\p^0 + CoAllullp^ej 
we finish the proof of Step 3. 
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Step 4: We remove the assumption that 6^ = c = 0 by moving the terms 
of bi and c to the right-hand side 


—ut -|- aijDijU — Xu = f — biDiU — cu. 

By the conclusion in Step 3 with = c = 0, there exists Aq = Xo{d, 5 ,p,9,Rq) 
such that for any A > Aq, 

+ 11 -^ ^Diu\\p^o 

< CoWf - biDiU - cuWp^e 

< Co\\f\\p,9 + CoKWDuWp^e + CoK\\u\\p^e. 


By taking A sufficiently large depending on d, 5, p, 6 , Rq, and K, we get 

ll^t||p,e + -^ll'*^llp,e + ^^ll-^^llp,6» + 11-^ ^ C'o||/||p, 0 - 

Step 5: To remove the assumption that T = oo, we simply follow the 
standard step in |Q, Theorem 2.1] or [1^, Theorem 6.4.1], Therefore, the 
estimate (^) is proved. 

Finally, in the case when bi = c = 0 and aij = aij{t), by using a scaling 
argument we can take Rq = 0. The theorem is proved. □ 
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